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An artificial two-atomic molecule, also called a 
double quantum dot (DQD), is an ideal system for 
exploring few electron physics [IlIllllllIlIIllSllllIIl 
[9l[10l[lTl[l2l[l3, 14J. Spin-entanglement between 
just two electrons can be explored in such sys- 
tems where singlet and triplet states are accessi- 
ble. These two spin-states can be regarded as the 
two states in a quantum two-state system, a so- 
called singlet-triplet qubit jl5| . A very attractive 
material for realizing spin based qubits is the car- 
bon nanotube (CNT) [16, 17, 18, 19, 20, 21j, be- 
cause it is expected to have a very long spin coher- 
ence time 53 [22l [23l [24l [25]. Here we show the 
existence of a gate-tunable singlet-triplet qubit 
in a CNT DQD. We show that the CNT DQD 
has clear shell structures of both four and eight 
electrons, with the singlet-triplet qubit present in 
the four-electron shells. We furthermore observe 
inelastic cotunneling via the singlet and triplet 
states, which we use to probe the splitting be- 
tween singlet and triplet, in good agreement with 
theory. 

Creating a qubit in a solid-state system demands con- 
trol of the number of interacting electrons. This control 
has to date been obtained using semiconducting materi- 
als operated close to the band-gap edge. We show in this 
Letter that shell structures in CNT DQDs, owing to the 
1-dimensional nature of the CNT, can be used to obtain 
the same kind of control. Both 4-electron and 8-electron 
(DQD)-shells are observed. We use one of the 4-electron 
shells to entangle the spin of two electrons and show that 
by separating the two electrons into separate QDs they 
form a spin triplet state, and by collecting them into 
the same QD they form a spin singlet state, i.e., a gate- 
tunable singlet-triplet qubit. 

The device analyzed in this Letter, schematically 
shown in Fig. [ija), is comprised of a CNT contacted 
by titanium electrodes, and gated by three top-gate elec- 
trodes, Gl, CG (center-gate), and G2, made of aluminum 
oxide and titanium. The device has two strongly coupled 
quantum dots in series as confirmed by the observation 
of the so-called honeycomb pattern in current (/gd) ver- 
sus voltage applied to Gl (Vbi) and G2 (Vg2) (Fig.[ljb) 
and Fig.[2|a)) [5 . The tunneling barrier between the two 
dots is due to a defect in the CNT (similar to Ref. [17 ). 
The resulting two dots have roughly equal charging ener- 
gies and level spacings (see below), from which we infer 
that the defect is located under or close to the center 



gate. The number of electrons in dot 1 and dot 2 can 
be controlled by tuning Vqi and Vq2^ respectively. In 
the middle of each hexagon (white areas in Fig. fl[b) and 
Fig.[2|a)) a fixed number of electrons are localized in each 
dot, and electron transport is suppressed by Coulomb 
blockade. Along the entire edge of the hexagons (blue 
lines), single electron transport is allowed through molec- 
ular states formed in the DQD, indicating a strong cou- 
pling between the two dots. The height (width) of the 
hexagons corresponds to the energy required to add an 
extra electron in dot 1 (dot 2). In Figjljb) the width and 
height of the hexagons alternate in size in a regular pat- 
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FIG. 1: Shell filling of a carbon nanotube double quan- 
tum dot. (a) Schematic figure of the device consisting of 
a carbon nanotube (CNT) contacted by titanium source and 
drain electrodes and gated with three top-gate electrodes, Gl, 
CG (center gate), and G2. Two quantum dots (QDl, QD2) 
are formed in series, one under Gl and one under G2. (b) 
Surface plot of current (/sd) at constant bias (V^d = 0.5 mV) 
as function of voltage applied to Gl (Vgi) and G2 (Vg2). 
Red numbers (N,M) are shell occupation numbers for one 
8-electron (DQD)-shell. Two further 8-electron shells were 
observed in connection the this shell, one below, and one to 
the right. 
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FIG. 2: Four-electron DQD-shell structure, (a) Surface 
plot of current (/sd) at constant bias (Fsd = 0.2 mV) as func- 
tion of voltage applied to Gl (Vgi) and G2 (Vg2)- The num- 
bers (N,M) are shell occupation numbers for one 4-electron 
shell, (b) Black lines: Schematic honeycomb diagram for a 4- 
electron shell with strong tunnel coupling, and a small cross 
capacitance. Gray lines: Same honeycomb diagram with neg- 
ligible tunnel coupling. Dashed lines indicate where the line- 
traces in Fig. |3] are measured, and the red oval indicates the 
region analyzed in Fig. ^ 



tern. The four hexagons marked with red numbers are 
distinctively larger than the other hexagons with three 
smaller hexagons in between, indicating that each dot 
has four-fold degenerate levels due to spin and orbital 
degeneracy [26l [27l [28j . An 8-electron shell structure of 
the DQD can therefore be identified in this plot. Shell 
occupation numbers (N,M), where N (M) is the level oc- 
cupation number in dot 1 (dot 2) are written onto the 
honeycomb diagram. 

The honeycomb diagram in Fig. [2|a) is measured for 
the same device but in another gate region where a new 
pattern in the sizes of the hexagons is observed. The 
hexagons alternate in size between large and small due 
to only spin degeneracy of the energy levels in each 
dot [26l [^, yielding a 4-electron shell structure of the 
DQD. The charging energies (/7ci, Uc2) and level spac- 
ings (A^i, AE2) for the two dots can be extracted 
from the honeycomb pattern as schematically shown in 
Fig.[2];b). The gate coupling of Gl (G2) to dot 1 (dot 2) 
is found from bias spectroscopy plots (not shown), and 
we find Uci ^ 3meV and AEi ^ 1.2 me V for dot 1, and 
Uc2 ^ 3.5 meV and AE2 ~ 1.5 meV for dot 2. Since 
charging energy and level spacing are almost identical 
for the two dots we deduce that the two dots are roughly 
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FIG. 3: Singlet and Triplet states in a four-electron 
shell, (a)-(b) Theoretical magnetic field dependence of the 
chemical potentials measured in (c)-(f). In (a) the 4-electron 
shell is filled from zero electrons (|00)), to one spin-down elec- 
tron in dot 2 (112)), to a singlet (Sb)- In (b) the device is 
filled from a spin-down electron in dot 1 (|ii)), to either a 
singlet (Sb) or a triplet (Tb-) (see text), to a three particle 
state with one spin-down electron in dot 1 and both a spin- 
down and spin-up electron in dot 2 (|iit2i2))- The line traces 
in (c)-(f) are extracted from honeycomb diagrams measured 
at B = 0,1,2,...,7T, where B is perpendicular to the nan- 
otube. (c) and (e) Horizontal and vertical line traces through 
hexagon (0,1) and (1,0) as function of electrostatic potential 
in dot 2 {E2) and dot 1 (^1), respectively, (d) and (f) Hor- 
izontal and vertical line traces through hexagon (1,1). Each 
line in (c)-(f) is offset 0.2 nA, 0.3 nA, 0.3 nA, and 0.5 nA re- 
spectively for clarity, and the left-most peak is centered at 



equal in size. We have observed both 4-electron and 8- 
electron shell structures in two different devices. 

We will in the rest of the Letter focus on a 4- 
electron shell with level spacings and charging energies 
similar the the 4-electron shell shown in Fig.[2|a), except 
AEi ~ 1.9 meV. The singlet ground state between region 
(1,1) and (0,2) is in general a bonding state of the local 
singlet (S(02), both electrons in dot 2), and the nonlocal 
singlet (S(ll), one electron in each dot): 



Sb =aS(ll)+/3S(02) 



(1) 



The detuning {e = E2 — Ei) dependent parameters a 
and [5 determine the weight of each state, and Ei and 
E2 are the electrostatic potentials in dot 1 and dot 2, 
respectively. Similarly for the triplets 



Tb- = a'T_(ll)+/3'T_(02) 
Tbo = a%(ll)+/3%(02) 
Tb+ = a'T+(ll)+/3'T+(02) 



(2) 
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FIG. 4: Singlet-triplet splitting probed by inelastic cotunneling. (a) Small section of the honeycomb diagram analyzed 
in Fig. p] with l^d = 50/iV at B = OT (left) and B = 6.5 T (right). The numbers (N,M) indicate electron occupation of the 
4-electron shell, (b) Chemical potentials for the singlet bonding (/xsb^oi) and triplet bonding (/xtb-^oi) with B = OT (solid 
green and blue lines) and with B = 6T (dashed green and blue lines) calculated using Eq. (|3| and Q. Onset of inelastic 
cotunneling which excites electrons from singlet to triplet (black arrow marked B) and from triplet to singlet (black arrow 
marked C) occurs when the separation between /iSe^oi and /xtb-^oi is equal to eV^d- (c) Schematic transport diagrams for 
elastic cotunneling (A) and inelastic cotunneling (B and C). (d) Surface plot of current (/sd) at l/gd = 0.2 mV versus magnetic 
field (B), and detuning (e) along the black dashed line in (a). Onset of inelastic cotunneling occurs along the white lines marked 
B and C, calculated using Eq. (51 with t = 0.32 meV, AE2 = 1.5 meV, V^a = 0.2 mV, and g — 2. Dashed grey lines indicate 
where the two inelastic cotunneling processes shown in (b) occurs. 



where — , 0, + denotes the spin magnetic moment in the z- 
direction, S^ = — 1,0,+1. We will in the following show 
the existence of the singlet and triplet states, i.e., the 
singlet-triplet qubit, on the basis of a magnetic field spec- 
troscopy on the 4-electron shell. 

In Fig. p[c) we analyze the magnetic field dependence 
of the width of hexagon (0,1), which involves 0, 1, and 2 
electrons in dot 2, The chemical potential for these two 
Coulomb peaks is given by p^: /ioi^oo oc —^qijl^B and 
/^Sb^oi oc -\-^g/j.BB, where /ioi^oo is the chemical po- 
tential for adding an electron to charge state (01) given 
no electron in the DQD-shell, and /iSe^oi is the chemi- 
cal potential for adding an electron in state Sb given one 
electron charge state (01). These two Coulomb peaks 
are therefore expected to separate by gfi^B as shown in 
Fig. Isja). The height of hexagon (1,0) are analogously 
expected to separate by gii^B. The measurements in 
Fig. [sjc) and (e) are in good quantitative agreement with 
the theory in Fig.jsja). The measured separation at 7 T is 
0.95 meV and 0.8 meV in (c) and (e), respectively, where 
theory predicts gfi^lT = 0.81 meV with ^ = 2 for nan- 
otubes. 

We now analyze the size of hexagon (1,1), which in- 
volves 1, 2, and 3 electrons in the DQD-shell. We show 
that by applying a magnetic field the 2-electron ground 
state can be changed from Sb to Tb- , which is used to 



estimate the exchange energy (J) (energy separation be- 
tween Sb and Tbo)- Transport at the first Coulomb peak 
in Fig.[3](d) is through different chemical potentials at low 
and high magnetic field, given by [13 //Sb^io ^ ^\9I^bB 
at low magnetic field {gfi^B < J), and /ixB-^io oc 
— \gi^BB at high magnetic field {gii^B > J). Simi- 
larly, transport at the second Coulomb peak in Fig. |3|d) 
is through /ii2^SB ^ —\gi^BB at low magnetic field 
{gji^B < J) and through /ii2^TB- ^ +^^Mb^ at high 
magnetic field {gji^B > J)[13 . The same magnetic 
field dependence is expected for the height of hexagon 
(1,1) (see Fig. |3jf)). Therefore, for increasing magnetic 
field, hexagon (1,1) decreases in size when Sb is ground 
state, and increases in size when Tb- is ground state, 
schematically shown in Fig. IsFb) . The measurements in 
Fig. IsFd) and (f ) are in good agreement with the theory 
in Fig.[3]^b) with the bend (shift of ground state from sin- 
glet to triplet) occurring at 5 ~ 2 - 3T, corresponding 
to an exchange energy of J ~ 0.23 - 0.35 meV. 

The exchange energy can, for large negative detun- 
ing (center of hexagon (1,1)), also be estimated from the 
tunnel coupling strength (t) using J c^ 4(tV^)^/^ci (see 
supplement material) [13l[30]. We estimate t ~ 0.32 meV 
from the curvature of the hexagons at the anticross- 
ings (see supplement material) [20 . This estimate of 
t yields a consistent estimate of the exchange energy 



J ^ A{t^/2f/Uci - 0.27 meV. 

The anticrossing between (1,1) and (0,2) (red area in 
fig. [2|b)) is analyzed in Fig. [ij We find that transport 
is governed by elastic and inelastic cotunneling via Sb 
and Tb- . The chemical potential for adding an electron 
to Sb and Tb- with Ei -\- E2 = 0, i.e., along the black 
dashed line in Fig. [4]^ a) is given by (see supplement): 
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We plot Eq. ^ and Q in Fig.[4];b) with B = OT (solid 
green and blue lines), and with B = 6T (dashed green 
and blue lines). We see that Sb is ground state for 

5 = OT, and that the two chemical potentials cross at 
elevated magnetic field, indicated with red arrow. 

At low magnetic field one broad peak in conductance 
versus detuning between (1,1) and (0,2) is seen (Fig.pFa), 
white arrow marked A). This conductance peak is due 
to elastic cotunneling via Sb, schematically shown in 
Fig. [4]^c) (mark A). Since elastic cotunneling via Sb in- 
volves both S(ll) and S(02), which have equal weight at 

6 = 0, the elastic cotunneling peak is centered around 
€ = 0. At high magnetic field the elastic cotunneling 
via Sb becomes suppressed because the ground state at 
e = changes from Sb to Tb- . Fig. Qd) shows a sur- 
face plot of /sd versus e and B along the black dashed 
line in Fig. Ilia). The white vertical line marked A is the 
expected position of the elastic cotunneling. 

At high magnetic field we observe two narrow peaks, 
marked B and C in Fig. [4Fa). These two narrow peaks 
are due to the onset of inelastic cotunneling via Sb and 
Tb- , schematically shown in Fig. llFc) mark B and C. 
Onset of inelastic cotunneling via Sb and Tb- occurs 
when the energy separation between their chemical po- 
tentials becomes equal to the applied bias: 
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We have from these two conditions calculated the onset 
of inelastic cotunneling in (e, 5)-space and plotted it as 
white lines marked B and C in Fi g. [4| ^d). Note that no 
fitting parameters are used in Fig. [4[d), the parameters 
used, t = 0.32 meV, AE2 = 1.5 meV were found in the 
analysis above. 



METHODS 

Fabrication and measurement setup 

The devices are made on a highly doped silicon sub- 
strate with a top layer of silicon dioxide. The CNTs are 
grown by chemical vapor deposition from islands of cat- 
alyst material and subsequently contacted by 50 nm Ti- 
tanium source and drain electrodes. Next, three narrow 



top-gate electrodes are fabricated between the source and 
drain electrodes, consisting of aluminum oxide and tita- 
nium [16] . A schematic figure of the device together with 
the measurement setup is shown in Fig. flla). Source- 
drain voltage (V^d) is applied to the source electrode and 
the drain electrode is connected through a current-to- 
voltage amplifier to ground. The three top-gate elec- 
trodes are named Gl, CG (center gate), and G2 starting 
from the source electrode. For the device reported on 
in this Letter we saw that Gl had a much lower gate- 
coupling than G2 and CG (see Fig. [ijb) and Fig. [2|a)) 
which we attribute to the Gl-electrode being damaged 
somewhere, weakening its gate-coupling. The gate cou- 
pling of Gl to dot 1 is aci = 2.9meV/V, and gate cou- 
pling of G2 to dot 2 is aG2 = 400 meV/V. The center gate 
is kept at FcG = OF for the measurements shown in this 
Letter. All data presented in this Letter are measured in 
a sorption pumped ^He cryostat at 350 mK. 
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